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The application of the theory of scale relativity to microphysics aims at recovering quantum 
mechanics as a new non-classical mechanics on a non-derivable space-time. This program was 
already achieved as regards the Schrodinger and Klein Gordon equations, which have been derived 
in terms of geodesic equations in this framework: namely they have been written according to a 
generalized equivalence/strong covariance principle in the form of free motion equations D 2 x/ds 2 = 
0, where D/ds are covariant derivatives built from the description of the fractal/non-derivable 
geometry Following the same line of thought and using the mathematical tool of Hamilton's bi- 
quaternions, we propose here a derivation of the Dirac equation also from a geodesic equation 
(while it is still merely postulated in standard quantum physics). The complex nature of the wave 
function in the Schrodinger and Klein-Gordon equations was deduced from the necessity to introduce, 
because of the non-derivability, a discrete symmetry breaking on the proper time differential element. 
By extension, the bi-quaternionic nature of the Dirac bi-spinors arises here from further discrete 
symmetry breakings on the space-time variables, which also proceed from non-derivability. 
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I. INTRODUCTION 



The theory of scale relativity generalizes to scale transformations the principle of relativity, which has been applied 
by Einstein to motion laws. It is based on the giving up of the assumption of space-time coordinate differentiability, 
which is usually retained as an implicit hypothesis in current physics. Even though this hypothesis can be considered 
£S) ' as roughly valid in the classical domain, it is clearly broken by the quantum mechanical behavior. It has indeed 
been pointed out by Feynman (see, e.g., Ref. [1]) that the typical paths of quantum mechanics are continuous but 
non-differentiablc. In the present paper, we apply the scale-relativistic approach to the understanding of the nature 
of the bi-spinors and of the Dirac equation in a space-time representation. 

One of the fundamental reasons for jumping to a scale-relativistic description of nature is to increase the generality 
i-^h ■ of the first principles retained. Einstein's principle of relativity relies on the assumption that coordinates are a priori 
differentiable. It is therefore not fully general, since it applies to coordinate transformations that are continuous and, 
at least two times, differentiable. The aim of the theory of scale relativity is to bring to light laws and structures that 
would be the manifestation of more general transformations, namely, continuous ones, either differentiable or not. 
The standard "general relativistic" theory will be automatically recovered as the special differentiable case, while new 
effects would be expected from the non-diffcrcntiable part. 



Giving up the assumption of differentiability has important physical consequences: one can show [2,3] that spaces 
of topological dimension Dt, which are continuous but non-differentiable, are characterized by a Dt measure which 
becomes explicitly dependent on the resolution (i.e., the observation scale) e at which it is considered and tends to 
infinity when the resolution interval e tends to zero. Therefore, a non-diffcrcntiable space-time continuum is neces- 
sarily fractal, in the general meaning given to this word by Mandelbrot [6], who thus means, not only self-similar, 
but explicitly scale-dependent in the more general way. This statement naturally leads to the proposal of a geometric 
tool adapted to construct a theory based on such premises, namely, a fractal space-time (see, e.g., Refs. [2,4,5,7-9]). 
Such an explicit and fundamental dependence of physical laws on the internal "resolution" scales leads to constrain 
the new scale laws (that need to be constructed) by a "principle of scale relativity" [2] . This principle extends to 
scale transformations (i.e., contractions and dilations of resolutions) the Galileo-Einstein principle of relativity, that 
was up to now applied only to displacements and motion. This allows to include the resolutions in the definition of 
the state of the reference system and to require scale covariance of the equations of physics under scale transforma- 
tions. A main consequence is that the geodesies of such a non-differentiable space-time are themselves fractal and in 
infinite number between any two points. This strongly suggests that the description could borrow some tools from 
the mathematical formalism of statistics, i.e., will provide a probabilistic interpretation, therefore naturally extending 
the realm of quantum behaviour to a larger spectrum of available scales. The domains of application are typically 
the asymptotic scale ranges of physics: (i) Microphysics: small length-scales and time-scales, (ii) Cosmology: large 
length-scales, (iii) Complex chaotic systems: long time-scales. 
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In earlier works [2,10] the Schrodinger equation and the (free) Klein-Gordon equation have been established, from 
the first principles of this theory, as geodesic equations on a fractal space/space-time. The third main evolution 
equation of standard quantum physics, the (free) Dirac equation is here derived, following the same trend of thought. 

As in these earlier works, the general method followed here, in order to reach a genuine demonstration of the Dirac 
equation from first principles, involves the following steps: (i) introduce a scale-space of internal resolutions, aimed 
at describing the underlying fractal structures of geodesic trajectories (that are identified with the "particles" in this 
approach) ; (ii) perform this description by writing scale differential equations coming under the principle of relativity, 
then solving them; (iii) express the effects induced on the laws of motion (space-time displacement) by the internal 
fractal structures of the underlying scale-space in terms of a covariant derivative ct/ds; when applied on the classi- 
cal position variable, this covariant derivative changes the nature of the classical velocity, which becomes complex 
(Schrodinger and Klein-Gordon equations) , quaternionic (Pauli equation), then bi-quatcrnionic (Dirac equation) at 
various depth levels of the theory; (iv) use the principle of relativity in the form of a generalized equivalence/strong 
covariance principle in order to write the motion equation as a geodesic equation, i.e. as a free motion equation of 
the formi^ x^/ds 2 = 0; (v) expand the covariant derivative and make a change of variable in order to introduce the 
quantum tool (the wave function is defined as a rc-cxpression of the action), then integrate and obtain the fundamental 
equations of quantum mechanics. 

In Sec. II, we actualize the former works dealing with this issue, proposing a more accurate and profound interpre- 
tation of the nature of the transition from the non-diffcrcntiable (fractal scales) to the differentiable (classical scales) 
domain, and carefully justifying the different key choices made at each main step of the reasoning. Then we derive, 
in Sec. Ill, the bi-quaternionic form of the Klein-Gordon equation, from which the Dirac equation naturally proceeds, 
as shown in Sec. IV. Section. V is devoted to the conclusion. 



II. THE FRACTAL TO CLASSICAL SCALE TRANSITION REVISITED 



A. Scale invariance and Galilean scale relativity 

A power law scale dependence, frequently encountered in natural systems, is described geometrically in terms of 
fractals [6,12], and algebrically in terms of the renormalization group [13,14]. As we show below, such simple scale- 
invariant laws can be identified with the "Galilean" version of scale-relativistic laws. 

Consider a non-differentiable (fractal) curvilinear coordinate £(x,e), that depends on some space-time variables 
x and on the resolution e. Such a coordinate generalizes to non-diffcrcntiable and fractal space-times the concept 
of curvilinear coordinates introduced for curved Riemannian space-times in Einstein's general relativity [2]. Rather 
than considering only the strict non-differentiable mathematical object C{x), we are interested in its different versions 
at various resolutions e. Such a point of view is particularly well-adapted to applications in physics, as any real 
measurement is always performed at a finite resolution. In this framework, C(x) becomes the limit when e — > of the 
family of functions C(x, e). But while C(x,0) is non-differentiable, C(x,e), which we call a fractal function [9,15], is 
differentiable for all e ^ 0. The physics of a given process will therefore be completely described if we can know £(x, e) 
for all values of e. C{x, e), being differentiable when e ^ 0, can be the solution to differential equations involving the 
derivatives of C with respect to both x and e. 

Let us now apply an infinitesimal dilation e — > e' = e(l + dp) to the resolution. Being, at this stage, interested in 
pure scale laws, we omit the x dependence in order to simplify the notation and obtain, at first order, 

C(e') = C(e + edp) = C(e) + ^j^edp = (1 + Ddp)C(e), (2.1) 

where D is, by definition, the dilation operator. The identification of the two last members of this equation yields 

This well-known form of the infinitesimal dilation operator shows that the "natural" variable for the resolution is 
lne, and that the expected new differential equations will involve quantities such as d£(x, e)/dlne. An example of 
equations describing scale dependence is given by the renormalization group equations, as first proposed by Wilson 
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[13,14]. 



We limit the present approach to the consideration of the simplest form that can be exhibited by this class of 
equations, leaving to future works the study of more complete cases. Let us built such an equation by stating that 
the variation of C under an infinitesimal scale transformation dine depends only on C itself; namely C determines 
the whole physical behavior, including the behavior under scale transformations. We thus write 



dC(x, e) 
9 In e 



/?(£). 



(2.3) 



As we are interested in the simplest form for such an equation, we expand /3(C) in powers of C. This can always 
be done since C may be renormalized, dividing it by its largest value, in such a way that the new variable C remains 
<C 1 in its variation domain. We obtain, to the first order, the linear equation 



of which the solution is 



dC(x, e) 
dine 



C(x, e) = Co(x) 



= a + bC , 



1 + t(x) 



(2.4) 



(2.5) 



where A b C( x ) is an integration constant and Cq = —a/b. 

For b < 0, Eq. (2.5) gives a fractal scale-invariant behavior at small scales (e <C A), with a fractal dimension 
D = D T — b. For a curvilinear coordinate C, the topological dimension D T is equal to 1. The scale dimension, 
S = D — Dt, is defined, following Mandelbrot [6,12], as 



dlnC 
~ dln(X/e) ' 

When S is constant, one obtains an asymptotic power law resolution dependence 



(2.6) 



C(x,e) = C (x) 



(2.7) 



Let us now check that for such a simple self-similar scaling law the new scale laws to be constructed come indeed 
under the principle of relativity extended to scale transformations of the resolutions e, as we have stressed in the 
introduction. The structure of the corresponding group of scale transformations can be established from the fact that 
the involved quantities transform, under a scale transformation e — > e', as 



ln ^l=ln^ 



Co 



Co 



(5(e) In ■ 



(2.8) 



6(e')=S(e). 



(2.9) 



These transformations have exactly the structure of the Galileo group, as confirmed by the dilation composition 
law, e — > e' — > e", which writes 



e" , e' , e" 
In — = In — + In — 

e e e' 



(2.10) 



It is worth noting that Eq. (2.5) gives also a transition from a fractal to a non-fractal behavior beyond some 
transition scale A. 
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B. Transition from non-differentiability (fractal scales) to differentiability (classical scales) 



The first consequence of the giving up of the coordinate differentiability is the differential proper time symmetry 
breaking. Relativistic motion involves a four-dimensional space-time with the proper time s as a curvilinear parameter. 
Strictly, the non-differentiability of the coordinates implies that the four- velocity 

v» = = hm + (2 . n) 

ds ds^o ds 

is undefined. This means that, when ds tends to zero, either the ratio dX^/ds tends to infinity, or it fluctuates 
without reaching any limit. However, as recalled in the introduction, continuity and non-diffcrcntiability imply an 
explicit scale dependence of the various physical quantities, and therefore of the velocity, V^. We thus replace the 
differential, ds, by a scale variable, 5s, and consider V fJ '(s,Ss) as an explicit fractal function of this variable. For a 
constant fractal dimension D, the resolution in X^, e^Ss), corresponds to the resolution in s, Ss, according to 

e M « Ss 1 / . (2.12) 

The advantage of this method is that, for any given value of the resolution, Ss, differentiability in s is recovered, 
which allows to use the differential calculus, even when dealing with non-differentiability. 



The scale dependence of the velocity suggests that we complete the standard equations of physics by new differential 
equations of scale. We therefore apply to the velocity and to the differential element, now interpreted as a resolution, 
the reasoning applied to the fractal function C in Sec. HA. Writing the simplest possible equation for the variation 
of the velocity V^(s, ds) in terms of the new scale variable ds, as a first order renormalization group-like differential 
equation of the form of Eq. (2.3); then, Taylor-expanding it, as in Eq. (2.4), using the fact that < 1 (in motion- 
relativistic units c = 1), we obtain the solution 



1 + <"(*) 



(2.13) 



where we have set b = 1/D — 1. Here, is the "classical part" of the velocity, Cl{V^), (see below the definition 
of the classical part operator C£), w M is the explicitly scale-dependent "fractal part" and r and £ M are chosen such 
that Ci(C,^) — and Ct((^) 2 ) = 1. We recognize here the combination of a typical fractal behavior, with a fractal 
dimension D, and of a breaking of the scale symmetry at the scale transition r, which is identified with the Compton 
scale of the system (r = h/mc), since w v^, when ds ^> t (classical behavior), and V* w w M , when ds <C r 
(fractal behavior). Recalling that D = 2 plays the role of a critical dimension, we stress that, in the asymptotic 
scaling domain, cx {ds/r)^ 1 / 2 [2], in agreement with Ref. [1], which allows to identify the fractal domain with the 
quantum one. In the present paper, only these simplest scale laws with a fractal dimension D = 2 are considered. 



The above description strictly applies to an individual fractal trajectory. Now, one of the geometric consequences 
of the non-differentiability and of the subsequent fractal character of space-time itself (not only of the trajectories) is 
that there is an infinity of fractal geodesies relating any couple of points of this fractal space-time. It has therefore 
been suggested [7] that the description of a quantum mechanical particle, including its property of wave-particle 
duality, could be reduced to the geometric properties of the set of fractal geodesies that corresponds to a given state 
of this "particle". In such an interpretation, we do not have to endow the "particle" with internal properties such 
as mass, spin or charge, since the "particle" is not identified with a point mass which would follow the geodesies, 
but its internal properties can simply be defined as geometric properties of the fractal geodesies themselves. As a 
consequence, any measurement is interpreted as a sorting out of the geodesies of which the properties correspond to 
the resolution scale of the measuring device (as an example, if the "particle" has been observed at a given position 
with a given resolution, this means that the geodesies which pass through this domain have been selected) [2,7]. 

The transition scale appearing in Eq. (2.13) yields two distinct behaviors of the system (particle) depending on the 
resolution at which it is considered. Equation (2.13) multiplied by ds gives the elementary displacement, dX^, of the 
system as a sum of two terms 

dX» = dx" + d£", (2.14) 
d£ M representing the "fractal part" and dx^ the "classical part" , each term being defined as 
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dx» = v" ds, 



(2.15) 



d^a^V^ds 2 ) 1 / 213 , (2.16) 

which becomes, for D = 2 

d£ M = a M V2Pds 1/2 , (2.17) 

with 21? = t. We note, from Eqs. (2.14) to (2.17), that dx^ scales as ds, while d£ M scales as ds? . Therefore, the 
behavior of the system is dominated by the d£ M term in the non-differentiable "fractal" domain (below the transition 
scale), and by the dx^ one in the diffcrcntiable "classical" domain (above the transition scale). 

Now, the Klein-Gordon and Dirac equations give results applying to measurements performed on quantum objects, 
but realised with classical devices, in the differentiable "classical" domain. The microphysical scale at which the 
physical systems under study are considered induces the sorting out of a bundle of geodesies, corresponding to the 
scale of the systems, while the measurement process implies a smoothing out of the geodesic bundle coupled to 
a transition from the non-differentiable "fractal" to the differentiable "classical" domain. We therefore define an 
operator C£( ), which we apply to the fractal variables or functions each time we are drawn to the "classical" domain 
where the ds behavior dominates. The effect of Ct is to extract, from the fractal variables or functions to which it is 
applied the "classical part", i.e., the part scaling as ds. 



C. Differential proper time symmetry breaking 

Another consequence of the non-differentiable nature of space-time is the breaking of the local differential proper 
time reflection invariance. The derivative with respect to the proper time s of a differentiable function / can be 
written twofold 

4f = lim m = lim M-fjs-ds) 

ds ds^O ds ds— >o ds 

The two definitions are equivalent in the differentiable case. One passes from one to the other by the transformation 
ds <-> — ds (local differential proper time reflection invariance), which is therefore an implicit discrete symmetry of 
differentiable physics. In the non-differentiable situation, both definitions fail, since the limits arc no longer defined. 
In the new framework of scale relativity, the physics is related to the behavior of the function during the "zoom" 
operation on the proper time resolution Ss, identified with the differential element ds. Two functions f' + and f'_ are 
therefore defined as explicit functions of s and ds 

f + (s, ds) = n" + ds,ds)-f{s,ds) (2 1Q) 



/' (s, ds) = n«M-f(s-ds,ds) 
ds 

When applied to the space-time coordinates, these definitions yield, in the non-differentiable domain, two four- 
velocities of which the components are fractal functions of the resolution, V+{x^(s),s,ds} and V^[x^(s), s, ds}. In 
order to go back to the "classical" domain and to derive the "classical" velocities appearing in Eq. (2.15), we smooth 
out each fractal geodesic in the bundles selected by the zooming process with balls of radius larger than r. This 
amounts to carry out a transition from the non-differentiable to the differentiable domain and leads to define two 
"classical" velocity fields now resolution independent: V_^[x >1 (s), s,ds > r] = Cl{V^[x fl (s), s, ds]) = v^[x M (s),s] and 
V^[x li (s),s,ds > t] = C£(V^[x fJ -(s), s, ds}) = [x^(s), s]. The important new fact appearing here is that, after the 
transition, there is no longer any reason for these two velocities to be equal. While, in standard mechanics, the 
concept of velocity was one- valued, we must introduce, for the case of a non-differentiable space, two velocities instead 
of one, even when going back to the "classical" domain. This two-valucdncss of the velocity vector finds its origin in 
the breaking of the discrete proper time reflection invariance symmetry (ds <-> — ds). Therefore, if one reverses the 
sign of the proper time differential clement, becomes v^_. But the choice of a reversed proper time element (—ds) 
must be as qualified as the initial choice (ds) for the description of the laws of nature, and we have no way, at this 
level of the description, to define what is the " past" and what is the "future" . The only solution to this problem is 
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to consider both the forward (ds > 0) and backward (ds < 0) processes on the same footing. Then the information 
needed to describe the system is doubled with respect to the standard diffcrcntiable description. 

A natural way of accounting for this doubling consists in using complex numbers and the complex product. This 
is the origin of the complex nature of the wave function of quantum mechanics, since this wave function can be 
identified with the exponential of the complex action that is naturally introduced in this framework [2,10]. The choice 
of complex numbers to represent the two-valuedness of the velocity is not an arbitrary choice. The use of these num- 
bers is the most natural way to generalize to two dimensions the set of real numbers, and obtain, in the relativistic 
motion case, the complex Klein-Gordon equation [10]. In the same way, as we shall see in Sees. Ill and IV, the use of 
bi-quaternionic numbers (i.e., complex Hamilton's quaternions) is the most natural way to generalize them to eight 
dimensions [16] and thus obtain the bi-quaternionic Klein-Gordon equation, which yields the Dirac equation. 



III. BI-QUATERNIONIC KLEIN-GORDON EQUATION 

It has long been known that the Dirac equation naturally proceeds from the Klein-Gordon equation when writen 
in a quaternionic form [17,18], i.e., using the quaternionic formalism, as introduced by Hamilton [19], and further 
developed by Conway [18,20]. However, this remains essentially formal in the standard framework, since there is no 
first principle reason for which the probability amplitude should be quaternionic. In the present approach, we shall see 
that its bi-quaternionic nature can be established as a consequence of the nonderivable nature of space-time. Indeed, 
we propose in the current section to introduce naturally a bi-quaternionic covariant derivative operator, leading to 
the definition of a bi-quaternionic velocity and wave- function, which we use to derive the Klein-Gordon equation in 
a bi-quaternionic form. This allows us to obtain, in Sec. IV, the Dirac equation as a mere consequence, since Dirac 
spinors and bi-quaternions are actually two equivalent representations of an electron. 



A. Further symmetry breakings and bi-quaternionic covariant derivative operator 

Most of the approach described in Sec. II remains applicable. However, the main new features obtained in the case 
we now study proceed from a deeper description of the scale formalism, considering the more general case when the 
peculiar choice of an identification of the differentials and the resolution variables (remember the choice explicited in 
Sec. II B, where we have set ds = Ss) is given up, implying the subsequent breaking of the symmetries: 

ds <-» — ds dx 11 <-> — dx 11 x^ <-» — x^ 




FIG. 1. Fractal function V M (s,e M ,e s ) 
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In the scaling domain, the four space-time coordinates X^(s, e M , e s ) are now four fractal functions of the proper time 
s and of the resolutions for the coordinates and e s for the proper time. We consider the case when, for an elemen- 
tary displacement dX^ corresponding to a shift ds in the curvilinear parameter, the resolutions verify e M < dX^ and 
e s < ds, which implies that, at a given X^, a forward shift ds of s yields a displacement dX M of X 11 and a backward 
shift — ds produces a displacement —dX^, of which the amplitudes are not necessarily equal (see Fig. 1). 

We can therefore apply, to these two different elementary displacements, the canonical decomposition 

dX" = dx" + df", with C£(dX") = dx» = vl ds, 

d£" = a^V^D{ds 2 )^, Cl(a^) = 0, Cl(«) 2 ) = 1, (3.1) 

and 

-dX^ = 5x» + S^, with a(-dX^) = Sx^ = v» ds, 

5? = a^V^Dids 2 )^ , Cl{at) = 0, Cl((a^) 2 ) = 1. (3.2) 

In the differcntiable case, dX^ = —(—dX^), and therefore v + = — vt . This is no longer the case in the non- 
diffcrentiable case, where the local symmetry dx^ <-> — dx^ is broken. 

Furthermore, we must also consider the breaking of the symmetry ds «-> — ds proceeding from the twofold definition 
of the derivative with respect to the curvilinear parameter s. Applied to X^ , considering an elementary displacement 
dX^, the classical part extraction process gives two classical forward and backward derivatives d/ds + and d/ds-, 
which yield in turn two classical velocities, which we denote v± + . Considering now the same extraction process, 

applied to an elementary displacement —dX' 1 , the classical forward and backward derivatives d/ds + and d/ds- allow 
again to define two classical velocities, denoted v± _ . We summarize this result as 

«t + =^, <-=^- (3-3) 
s f! ds± s n ds± 

Contrary to what happens in the diffcrentiable case, the total derivative with respect to the proper time of a fractal 
function f(x(s),s) with integer fractal dimension contains finite terms up to higher order [21] 

df__d£ df_dX^_ 1 d 2 f dX^dX v 1 d 3 f dX x dX»dX v 
ds ds dx^ ds 2 dx^dx v ds 6 dx x dx^dx u ds 

In our case, a finite contribution only proceeds from terms of D-order, while lesser-order terms yield an infinite 
contribution and higher-order ones are negligible. Therefore, in the special case of a fractal dimension D = 2, the 
total derivative writes 

<£_d£ dX» 1 d 2 .f dX^dX v 

ds~ ds + ^' ds + 2 dx^dx" ds ( ' ' 

Usually the term dX^dX v jds is infinitesimal, but here its classical part reduces to C£(d^ d£")/ds. 

We can, at this stage, define several total derivatives with respect to s of a fractal function /. We write them - using 
a compact straightforward notation with summation over repeated indices - after substituting, in the four-dimensional 
analog of Eq. (3.5), the expressions for the derivatives of the X' 1 obtained when using Eqs. (3.1) and (3.2) with the 
different expressions of Eq. (3.3) for the classical velocities 

^ ± ± ± ± ± = ^ + « ± + < ± + a^alV-^l— , (3.6) 

ds s * v z t ds s i± s k aX^ dX^dX v 

with w» = a^V2D{ds 2 )^-^. (3.7) 

Now, when we apply the classical part operator to Eq. (3.6), using Eq. (3.7) and the properties of the a± 's as 
stated in Eqs. (3.1) and (3.2), the 's disappear at the first order, but, at the second order, for the fractal dimension 
D = 2, the fractal behavior of Eq. (2.16) writes 
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C£(w1 ± w" ±± ) = TZVv^ds, (3.8) 

the =F sign in the right-hand side being the inverse of the s-sign in the left-hand side. 

Equation (3.8) follows from the fact that the c?£ M 's are of null classical part and mutually independent. The 
Minkowski metric component, rf v implies indeed that the classical part of every crossed product w± ±u>± ±, with 

S /J, S V 

\i 7^ v, is null. This is due to the fact that, even if each term in the product scales as ds 1 / 2 , each behaves as an 
independent fractal fluctuation around its own classical part. Therefore, when we proceed to the smoothing out 
of the geodesic bundle during the transition from the fractal to the classical domain, we apply a process which is 
mathematically (not physically) equivalent to a stochastic "Wiener" process, and also more general, since we do not 
need any Gaussian distribution assumption. Thus, we can apply to the classical part of the w crossed product the 
property of the product of two independent stochastic variables: i.e., the classical part of the product is the product 
of the classical parts, namely here zero. 

Thanks to Eq. (3.8), the last term of the classical part of Eq. (3.6) amounts to a Dalambertian, and we can write 

f s1t ,^ = (§- s+ v^ t d^V d ^)f, (3.9) 

where the =F sign in the right-hand side is still the inverse of the s-sign. When we apply these derivatives to the 
position vector X* 1 , we obtain, as expected, 

^ ±± =< ± . (3.10) 

US » c s f< 

We consider now the four fractal functions — X^(s, e^, e s ). At this description level, there is no reason for 
(— X^)(s, e M , e s ) to be everywhere equal to — (X^)(s, e M , e s ), owing to a local breaking of the P (for \i = x,y,z) 
and T (for /1 = t) symmetries. Furthermore, as we have stressed above for X^, at a given —X^, a forward shift ds of 
s yields a displacement d(-X^) of -X' 1 and a backward shift — ds produces a displacement — d(— X^), with no nec- 
essarily equal amplitudes. We can therefore apply, to these two different elementary displacements, a decomposition 
similar to the one described in Eqs. (3.1) and (3.2), i.e., 

d(-X^) = dx" + df", with Ci{d{-X»)) = dx» = vl ds, 

d> = a^V2p(ds 2 )*, ««} = o, a((d t i) 2 ) = i, (3.11) 

and 

-d(-X") = fa" + tff, with Ci{-d{-X»)) = fix" = vt ds, 

5^ = ~atVZD{ds 2 )^ j C£(a^)=0, «((5^) 2 ) = 1, (3.12) 

Then, we jump to the classical parts and we are once more confronted to the breaking of the ds <-> —ds symmetry. 
Applied to — A", considering an elementary displacement d(— A"), the classical forward and backward derivatives 
d/ds+ and d/ds- give again two "classical" velocities, which we denote t;± + . After another classical part extraction, 

S fl 

considering now an elementary displacement — d(— A"), the classical forward and backward derivatives d/ds+ and 
d/ds- yield two other "classical" velocities, denoted v± _ . We write in short 

s 

~n dx^ i ■ Sx^ . . 

vl + =— , 5£_=_. (3.13) 

s fi as± s f. ds± 



We therefore obtain new different total derivatives with respect to s of a fractal function /, which we write 

T- ± ± ± ± ± = if + ± + < ± )4k + '<~ a ± v 7^f^ > ( 3 - 14 ) 
ds s x v * t ds s f. s ^ oX^ dX^dX v 

with tip = d^VZDids 2 )^-^. (3.15) 
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The classical part operator applied to Eq. (3.14) yields the classical total derivatives 



df 
ds 



± ± ± ± ± 



d_ 

d~ s 



When wc apply these derivatives to the position vector X 11 , we obtain, as expected again, 

dX" 



ds 



(3.16) 



(3.17) 



If we consider the simplest peculiar case when the breaking of the symmetry dx^ <-> —dx^ is isotropic as regards the 
four space-time coordinates (i.e., the signs corresponding to the four fi indices are chosen equal), we are left with eight 
non-degenerate components - four v± ± and four v±±- which can be used to define a bi-quaternionic four-velocity, 



v" = i« + +^_)-i« + -(£_) + 



«--<+) 



ei 



+ 



{vt_+vl_)-Uvt_-vt + ) 



e 2 + 



±(i£ + + t£ + )-i(t£ + + t5£_) 



e 3 - 



X ~~ 2 

with a(i — 1,2,3) denoting Hamilton's imaginary units satisfying the associative but non-commutative algebra 



(3.18) 



GiEj — 5ij + ^ ' £jjfe6fc, i,j — 1, 2, 3 



(3.19) 



k=l 



where e^fc is the usual completely antisymmetric three- index tensor with £123 = 1. 



The freedom in the choice of the actual expression for V M will be discussed later. It is constrained fy the fol- 
lowing requirements: at the limit when e M — ► dX^ and e s — > ds, every e^-term in Eq. (3.18) goes to zero, and, 
as = vt + in this limit, one recovers the complex velocity of the scale-relativistic Schrodinger equation [2], 
V = [v+ + + v' J i + — i(v+ + — vt + )\/2; at the classical limit, every term in this equation vanishes, save the real term, 
and the velocity becomes classical in its usual meaning, i.e., real: 



J ++- 



The bi-quaternionic velocity thus defined corresponds to a bi-quaternionic derivative operator d /ds, similarly de- 
fined, and yielding, when applied to the position vector X^, the corresponding velocity. For instance, the derivative 
operator attached to the velocity in Eq. (3.18) writes 



+ 



i. 

ds 

1 

2 



1 

" 2 

' d_ 
ds ' 



d_ d_ 

ds~ +++ ds~~ 



ds 



1 

' 2 

' d_ 

ds' 



d_ 

ds 



++ 



ds ' 



d_ 

ds ' 



+ 



d_ 

ds 

d_ 

ds 



d_ 

ds ' 



d_ 

ds 



d_ 

ds 



++ 



ds 



++ 



ild d 
2 1 ds~~ + + 



ei 



e3- 



(3.20) 



Substituting Eqs. (3.9) and (3.16) into Eq. (3.20), we obtain the expression for the bi-quaternionic proper-time 
derivative operator 



[l + ei + e 2 + (l 



d 



(3.21) 



the + sign in front of the Dalambertian proceeding from the choice of the metric signature (+,—,—,—). We keep 
here, for generality, the d/ds term, stressing that it is actually of no use, since the various physical functions are not 
explicitly depending on s. It is easy to check that this operator, applied to the position vector X^, gives back the 
bi-quaternionic velocity V M of Eq. (3.18). 

It is worth noting that the expression we have written for in Eq. (3.18) is one among the various choices we 
could have retained to define the bi-quaternionic velocity. The main constraint limiting this choice is the recovery 
of the complex and real velocities at the non-relativistic motion and classical limits. We also choose V M such as 
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to obtain the third term in the right-hand side of Eq. (3.21) under the form of a purely imaginary Dalembertian, 
which allows to recover an integrable equation of motion. To any bi-quaternionic velocity satisfying both prescriptions 
corresponds a bi-quaternionic derivative operator if /ds, similarly defined, and yielding this velocity when applied to 
the position vector X^. But, whatever the definition retained, the derivative operator keeps the same form in terms 
of the bi-quaternionic velocity V M , as given by Eq. (3.27). Therefore, the different choices allowed for its definition 
merely correspond to different mathematical representations leading to the same physical result. 



B. Bi-quaternionic wave-function 

Since V M is bi-quaternionic, the lagrange function is also bi-quaternionic and, therefore, the same is true of the 
action. 

A generalized equivalence principle, as well as a strong covariance principle, allows us to write the equation of 
motion under a free-motion form, i.e., the differential geodesic equation 

^ = 0, (3.22) 

ds 

where is the bi-quaternionic four-velocity, e.g., the covariant counterpart of V* defined in Eq. (3.18). 
The elementary variation of the action, considered as a functional of the coordinates, keeps the usual form 

5S = -mcVpSx* 1 . (3.23) 
We thus obtain the bi-quaternionic four-momentum, as 

Vp = mcV^ = -dpS. (3.24) 
We are now able to introduce the wave function. We define it as a re-expression of the bi-quaternionic action by 

1>-%1> = ^r^«S, ( 3 - 25 ) 
using, in the left-hand side, the quaternionic product derived from Eq. (3.24), as 

y = j^V" 1 ^- (3.26) 

m 

This is the bi-quaternionic generalization of the definition used in the Schrodinger case: ip = e lS / So . It is worth 
stressing here that we could choose, for the definition of the wave function in Eq. (3.25), a commutated expression in the 
left-hand side, i.e., (d^ip)^ 1 instead of d^ip. But with this reversed choice, owing to the non-commutativity of the 
quaternionic product, we could not obtain the motion equation as a vanishing four-gradient, as in Eq. (3.33). Therefore, 
we retain the above simplest choice, as yielding an equation which can be integrated. This non-commutativity induced 
property might be considered as a key for a future understanding of the parity violation, which will not be developed 
here. 



C. Free-particle Klein-Gordon equation 

As, in what follows, we only consider s-stationary functions, i.e., functions which do not explicitly depend on the 
proper time s, the derivative operator reduces to 

^- = Vd v + iVd u d v . (3.27) 
ds 

Now this expression is independent of the peculiar choice retained for the bi-quaternionic form of the four-velocity, 
as the representation dependent term in the right-hand side of Eq. (3.21) has vanished. 

The equation of motion, Eq. (3.22), thus writes 
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iy v d v + ivd v d„) v p = o. (3.28) 

We replace V^, respectively V v , by their expressions given in Eq. (3.26) and obtain 

i^l ( i^~ l d v i)d v + iVd v d^\ U^d.if) = 0. (3.29) 
m \ m J v ' 

The particular choice So = 2mD 1 allows us to simplify this equation and we get 

i>- l d v i> d v (if>-%^) + ^8^-%^) = 0. (3.30) 
The definition of the inverse of a quaternion 

^V _1 =^"V = 1, (3-31) 

implies that tp and ip^ 1 commute. But this is not necessarily the case for ip and d^ip -1 nor for ip" 1 and d^ip and 

their contravariant counterparts . However, when we derive Eq. (3.31) with respect to the coordinates, we obtain 

w 1 = -(^v# _1 

1>-%il> = -(d^- 1 )^, (3.32) 
and identical formulae for the contravariant analogues. 

Developing Eq. (3.30), using Eqs. (3.32) and the property d v d v d il — d^d v d v , we obtain, after some calculations, 

duftVd^W- 1 ] = 0. (3.33) 

We integrate this four-gradient as 

(ff'd^ip- 1 + C = 0, (3.34) 
of which we take the right product by ip to obtain 

d v d u ip + Cip = 0. (3.35) 
We therefore recognize the Klein-Gordon equation for a free particle with a mass m, so that m 2 c 2 /h 2 = C. 



IV. DIRAC EQUATION 

We now use a long-known property of the quaternionic formalism, which allows to obtain the Dirac equation for a 
free particle as a mere square root of the Klein-Gordon operator (see, e.g., Refs. [17,18]). 

We first develop the Klein-Gordon equation as 

1 d 2 tp d 2 ip d 2 ip d 2 ip m 2 c 2 



dt 2 dx 2 ' dy 2 dz 2 h 2 



i>. (4.1) 



Thanks to the property of the quaternionic and complex imaginary units e\ = e 2 = e 2 = i 2 = — 1, we can write 
Eq. (4.1) under the form 



: At first sight, the meaning of the choice So = 2mT> is not a necessary condition from the viewpoint of physics, but merely a 
simplifying choice as regards the equation form. Indeed, it is only under this particular choice that the fundamental equation 
of dynamics can be integrated. However, if we do not make this choice, the ip function is a solution of a third order, non linear, 
complicated equation such that no precise physical meaning can be given to it. We therefore claim that our choice <So = 2mX> 
has a profound physical significance, since, in our construction, the meaning of ip is directly related to the fact that it is a 
solution of the Klein-Gordon and Dirac equations obtained below. 
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1 d 2 1p 2 d 2 1p 2 . 2 2 V> . , 2 , . 2 ^ 2 C 2 2,2 ,„ x 

? a* = e ^ 62 + ie >w l + e ^ ei + ' — e ^- (42) 

We now take advantage of the anticommutative property of the quaternionic units (e,ej = —ejei for i 7^ j) to add 
to the right-hand side of Eq. (4.2) six vanishing couples of terms which we rearrange as 

13/1 dtp\ d f dip dip . dtp mc , \ 

c dt {cat) = 63 to I 63 to 62 + 61 V + 63 to 61 - *X e3 ^ 3 J 62 

d f dip dip . dip mc 

+ ei-r- e 3 — e 2 + ei — z + e 3 — e x - i—e 3 ipe 3 
dy \ ox dy dz n 

d / dip dtp . dip mc , \ 

+ e 3 — e 3 — e 2 + e 1 — 1 + e 3 — e 1 - i—e 3 ipe 3 e x 
oz \ ox dy dz n ) 

mc ( dip dip . dip mc \ 

We see that Eq. (4.3) is obtained by applying twice to the bi-quaternionic wavefunction ip the operator 

19 d d d mc 

~"H7 = e 3^~ e 2 + e l^~ Z + e 3^- e l - l -T e 3\ ) e 3- (4-4) 

c ot to dy dz n 

The three first Conway matrices e 3 ( )e2, ei( )i and e 3 ( )ei [22], figuring in the right-hand side of Eq. (4.4), can 
be written in the compact form —a k , with 



<7fc 

a k 




the Ofc's being the three Pauli matrices, while the fourth Conway matrix 

e 3 ( )e 3 = 



can be recognized as the Dirac (3 matrix. We can therefore write Eq. (4.4) as the non-covariant Dirac equation for a 
free fermion 

19V k mC a 1 ,Ar\ 

^ = -^-'T W - (4 ' 5) 
The covariant form, in the Dirac representation, can be recovered by applying ie 3 ( )e 3 to Eq. (4.5). 

The isomorphism which can be established between the quaternionic and spinorial algebrae through the multipli- 
cation rules applying to the Pauli spin matrices allows us to identify the wave function ip to a Dirac spinor. Spinors 
and quaternions are both a representation of the SL(2,C) group. See Ref. [23] for a detailed discussion of the spinorial 
properties of bi-quaternions. 



V. CONCLUSION 



The last of the three most fundamental motion equations (Schrodinger, Klein-Gordon and Dirac) which are merely 
postulated in standard quantum mechanics has been established here, in the framework of Galilean scale relativity, 
as a geodesic equation in a fractal space-time. The change from classical to quantum relativistic motion arises from 
successive symmetry breakings in the fractal geodesic picture. 

First, the complex nature of the wave function is the result of the differential (proper) time symmetry breaking, 
which is the simplest effect arising from the fractal structure of space (space-time). At this stage, Galilean scale 
relativity with a complex wave function permits the derivation of both the Schrodinger [2] and Klein-Gordon [10] 
equations. 
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To go on with the description of the elementary properties encountered in the microphysical world, we have con- 
sidered here further breakings of fundamental discrete symmetries, due to nonderivability, namely the differential 
coordinate symmetry (dx^ <-> —dx 1 *) breaking and the parity and time reversal symmetry breakings. These new 
breakings provide a four-complex component wave function (i.e., a eight component wave function), of which the 
most natural mathematical representation is in term of bi-quaternionic numbers [16]. We therefore obtain the spino- 
rial and the particle anti-particle nature of elementary objects which we can describe as Dirac spinors. Here, spin 
arises from the isomorphism between the quaternionic and spinorial representations, both of which are different repre- 
sentations of the SL(2,C) group. At this stage, the Klein-Gordon equation is written in a bi-quaternionic form which 
naturally yields the free Dirac equation. 

It is worth stressing that these results only proceed from a restricted use of the scale-relativistic potentialities. We 
have, in the present work, limited our investigations to the induced effects of scale laws on the equations of motion, in 
the framework of dilation laws exhibiting a Galilean group structure, i.e., a fractal space-time with a constant fractal 
dimension D = 2. This is only one of the simplest levels at which the scale-relativistic program can be achieved. 
We have also made a series of simplifying choices, which we have explained and justified all along the derivation 
procedure. Some of them have been dictated by physical or experimental considerations, but others only correspond 
to special cases, provisionally retained, so that other possibilities will have to be explored in the future. 
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